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•^r ' A random Fibonacci sequence is defined by the relation f/„ = \gn-i ± <?n-2|, where 

the ± sign is chosen by tossing a balanced coin for each n. We generalize these 

sequences to the case when the coin is unbalanced (denoting by p the probability of a 

+), and the recurrence relation is of the form g„ = |A(7„_i ±gn-2\- When A > 2 and 

< p < 1, we prove that the expected value of g„ grows exponentially fast. When 

Pm i A = Afe = 2 cos(7r/fc) for some fixed integer fc > 3, we show that the expected value 

P^ ' of Qn grows exponentially fast for p > (2 — Afc)/4 and give an algebraic expression 

r^ . for the growth rate. The involved methods extend (and correct) those introduced 

td : in 0. 
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^ ', 1. Introduction 

^—v ■ A random Fibonacci sequence is a sequence {gn)n defined by its first two terms gi 

QQ I and 52 (which in the sequel arc assumed to be positive) and the recurrence relation 

^^ ■ 5n+i — \gn ± .9n-i|, whcrc for each n the ± sign is chosen by tossing a balanced 

coin. 
- A generalization of this notion consists in choosing the ± sign by an unbalanced 

/\ ' coin (say: + with probability p and — with probability q := 1—p). Another possible 

^4 • generalization consists in fixing two real numbers, A and /^, and considering the 

recurrence relation g„ = |Ag„_i ± /ig„_2|, where the ± sign is chosen by tossing a 
balanced (or unbalanced) coin for each n. By considering the modified sequence 
5n '■= ^h/a*" , which satisfies 5„ ~ |-j=^„_i ± gn-2\: we can always reduce to the 
case fj. = 1. In the following, we refer to the random sequences g„ = \Xgn-i ±5„_2| 
where the + sign is chosen with probability p and the — sign with probability 
q := 1 — p a.s {p, X)-random Fibonacci sequences. 

In [3] , we investigated the question of the asymptotic growth rate of almost all 
{p, l)-random Fibonacci sequences. We obtained an expression of this limit which is 
simpler than the one given by Divakar Viswanath in j^ and which is not restricted to 
p = 1/2: It is given by the integral of the natural logarithm over a specific measure 
(depending on p) defined on Stern-Brocot intervals. The corresponding results for 
(p, A)-random Fibonacci sequences, which involve some techniques presented here 
but also complementary considerations, is the object of another publication (see [2]). 
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Here, we are concerned with the evaluation of the Uniit value of (E((7„))^/", where 
(<?n)n is a (p, A)-randoni Fibonacci sequence and E stands for the expectation. 

This limit was studied in the particular case p = 1/2 and A = 1 in [S], where 
the growth rate of the expected value of a (1/2, l)-random Fibonacci sequence is 
proved to be asymptotically equal to a — 1 ~ 1.20556943, where a is the only real 
zero of a'^ = 2a^ + 1. The proof involves the study of the binary tree T naturally 
defined by the set of all (1/2, l)-random Fibonacci sequences (if a node a has 6 as a 
child, then the node h has two children, labelled by a + 6 and by |a — 5|). The study 
of T is made by considering the biggest subtree of T (denoted by R) which shows 
no redundance, that is in which we never see two different edges with the same 
values a and h (in this order) as parent and child. This restricted tree R has many 
combinatorial and number-theoretic aspects which are of interest. Let us mention 
that, after the publication of [5], we realized that there was a problem in the last 
step of the proof of its main result; this mistake is explained in Remark [1] 

In the present paper, we correct and extend the result of [5j to some (p, A)- 
random Fibonacci sequences: For any p S [0, 1], any A > 2, and any A of the form 
2cos(7r/A:) (denoted by Afc), where A: > 3 is an integer. 

For A = Afe , the combinatorial properties of the tree R extends in a surprinsingly 
elegant way, leading to an extremely natural generalization of the results previously 
mentioned. In particular, the link made in [3] and [S| between the tree R and 
continued fraction expansion remains true for A^ = 2 cos(7r/fc) and corresponds to 
so-called Rosen continued fractions, a notion introduced by David Rosen in [6]. 
We will not resort to continued fractions in the present work, but this aspect is 
presented in 0. 

An interesting fact to notice is that the values A^ and A > 2 arc the only ones 
for which the group of transformations of the hyperbolic half plane H^ generated 
by the transformations z i — > ^l/-z and z i — > z -I- A, said to be a Hecke group, is 
discrete (cf. [l]). We will not use that fact in the following, but it suggests that it is 
highly probable that some link is to be made between random Fibonacci sequences 
and hyperbolic geometry; in particular, possible future extensions of the combina- 
torial point of view given by the restricted trees defined below could have some 
interpretation in hyperbolic geometry for values of A for which the corresponding 
Mobius group is not discrete. 

We are grateful to Kevin Hare, of University of Waterloo (Canada), for helpful 
comments, and to Jean-Frangois Quint, of CNRS and Universite Paris-13 (France), 
for fruitful discussions on hyperbolic geometry. 

2. Results 

Our main results are the following theorems. 

Theorem 1. Let p £ [0, 1], k > 3 and irin be the expected value of the n-th term 
of a {p, Xk) -random Fibonacci sequence. 

• If p > Pc := (2 - Afe)/4, then 

' ^^yp) 1 + — rv^ > 1' 

m„ «^oo \^ ak[pr J 

where a^ {p) is the only positive root of the polynomial 
Pk{X) := X^^ - XkX"^^-^ - (2p - l)X^^-'^ - \kpq^-^X^-^ - p^q^^-"^. 
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• If p = pc, then (m„)„ grows at most linearly. 

• If p < Pc, then {mn)n is bounded. 



l/n 

5" 



Hence the critical 



Note that by Jensen's inequality, we have E[g„]^/" > E 
value for the growth rate of the expected value is smaller than the critical value 
when one considers the almost-sure growth rate. It is proved in [2] that the latter 
is equal to 1/fc, which is strictly larger than pc = {2 — Afc)/4 = (1 — cos(7r/fc))/2. 

Theorem 2. Let A > 2, < p < 1, and rrin he the expected value of the n-th term 
of a {p, X)-random Fibonacci sequence. Then 



m„+i ^ A + y/A^ + 4(2p - 1) 

77l„ n^oo 2 

In view of the study of {p, A)-random Fibonacci sequences for other values of A, 
we also investigate an aspect of the regularity of the behaviour of the growth rate 
of such a sequence in the neighbourhood of A = 2. 

Corollary 1. Let < p < 1. If we assume that, for any A in the neighbourhood of 
2, the expected value of a {p, X) -random Fibonacci sequence increases exponentially 
fast with growth rate equal to G{X), then Q cannot be analytic at X = 2. 

In the particular case p = 1/2, we can even prove the non-analyticity of the 
growth rate on any left neighbourhood of 2. 

Corollary 2. Let p ~ 1/2. Assume that, for any X < 2, the expected value of a 
(\ / 2, X) -random Fibonacci sequence increases exponentially fast, with growth rate 
equal to Q{X). If Q is differentiable at X = 2, then G'{2) = 1. If Q is of class C" at 
A = 2, then t/(*)(2) = for any i e [2, n]. As a corollary, Q cannot be analytic at 
A = 2. 

We first prove Theorem [2] in Section [3l since the case A > 2 is much easier. 
Section|3]is devoted to general facts about the reduced tree for A = A^. It introduces 
notations and useful tools for the proof of Theorem [1] This theorem is proved in 
Sections [5l [6] and [71 by extending the ideas introduced in [5j for the case k = 
3 and p = 1/2. Proofs of corollaries [T] and [2] about the non-analyticity in the 
neighbourhood of 2 can be found in Section [H Section [9] contains a discussion 
about open questions. 

3. Case A > 2 

We introduce the tree T\{a,b), which shows all different possible {p, A)-random 
Fibonacci sequences with first positive terms gi ^ a and 92 = b: The root of the 
tree Tx{a, b) is labelled by a, its unique child is labelled by 6, and each node of the 
tree labelled by (3 and with parent labelled by a has exactly two children, the right 
one labelled by A/3 -I- a and the left one by |A/3 — a|. When (3 is the label of a node 
in Tx{a,b) with parent labelled by a, it will be convenient to consider the vector 
(a, (3) as the label of the corresponding edge. We also talk about the right and 
left children of an edge labelled (a, /3), which are respectively the edges labelled by 
(/3, A/3 -Ha) and (/3, |A/3 - a|). 

We also introduce the weight of an edge in the tree T\{a,b), corresponding to 
the probability that a (p, A)-random Fibonacci sequence passes through this edge. 
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^2 
^3 



11 111 31 5 

Figure 1. The beginning of the tree Ti(l, 1). 

More formally, the initial edge in T\{a,b) has weight 1. If an edge has weight w, 
then its right and left children have respective weight pw and qw. (Recall that 
q^l-p.) 

We organize the edges of the tree T\{a, b) in rows: The initial edge is the only 
edge in row 2; any child of an edge in row n is in row n + I. For any n > 2, we 
denote by ipn the set of edges in row n in T\{a, b). 

The average of any subset X of edges of T\{a, b) is defined as the sum M{X) of 
all terms of the form [3w, where /3 is the second coordinate of an edge in X and w 
is the weight of this edge. Observe that the expected value rUn of the n-th term of 
a (p, A)-random Fibonacci is given by M{ip„)- 

It will be of interest to consider the sequence {£s) of numbers read along the 
leftmost branch of T\{a, b): ii := a, £2 := b and €5+1 := l^^s — £s-i\ for s > 2. 

We now turn to the proof of Theorem [21 

Lemma 1. If b > a and X> 2, then for any edge in T\{a, b) labelled by {a, (3), we 
have P > a. 

Proof. This is immediately proved by induction. D 

Corollary 3. Ifb>aandX> 2, then M{ip„) = AA/('0„-i) + (2p - l)A/('0„-2)- 

Proof. We deduce from Lemma [T] that absolute values are never used in the com- 
putations of labels in Txio-, b). Regrouping the contributions of the two children of 
an edge in row n — 1, and summing over row 71 — 1, we get the result. D 

The preceding corollary shows that for b > a and A > 2, there exists constants 
C and C" (depending on a and b) such that 

M(V'„+2) = Ca" + C'a'", 
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where a := (A + ^X^ + A{2p - \))/2 and a' := (A - ^A2 + 4(2p - l))/2 arc the 
roots of the polynomial X^ — XX — {2p — 1). Observe that, by Lemma[Tl if 6 > a 
then all labels in the tree T\{a, h) are larger than &, hence Af (V;„+2) > 6 > 0. Since 
a > 1 and \a'\ < 1, we deduce that C > 0, and Theorem [2] is proved when b > a. 

li b < a, wc first consider the case where the sequence (Ig) along the leftmost 
branch is unbounded. Then there exists a first S > 2 such that Is+i > f-s- We 
can then consider T\{a,b) as the disjoint union of the trees T\{£s,X£s + ^s-i), 
2 < s < S, and the tree T\{is,^s+i)- M{^n+2) can be written as a convex 
combination of similar expressions in these trees. Since for each one of these trees, 
the labels of the first edge are well-ordered. Theorem [2] is valid for them. A simple 
computation shows that the result extends to T\{a,b). 

It remains to consider the case where (is) is bounded. We then consider T\{a, b) 
as the disjoint union of the leftmost branch and infinitely many trees, namely the 
trees Tx{is,Xis + ^s-i), 2 < s, whose first-edge labels are well-ordered. For each 
s, there exists constants Cs and C[. (depending on £s_i and is) such that 



n-l 
s=Q 



M(V„+2) = q"in+2 + V q'p (C,a"-^+i + Cy- 



Using the fact that Cs and C's arc bounded, we obtain that M(V'n+2) ^ Ka" for 
some K > 0, which ends the proof of Theorem [2] 

4. Reduced tree in the case A = Afe for some fc > 3 

From now on, we fix an integer fc > 3, and set A := Afc. We keep the notations 
concerning the tree TA(a,6) introduced in the preceding section. 
We introduce the two matrices 

a, z.(J -) a... ..(; I 

Observe that the right child of an edge labelled (a,/3) is labelled by {a,(3)R, and 
the left child is labelled by {u, \v\) where (w, v) = {a, (3)L. 

Definition 1. Any right child in a tree is said to be a 0-tli left child. For any 
integer m > 0, a child in a tree is said to be an m-th left child ijf it is the left child 
of an (m — l)-th left child. 

Proposition 1. Let {a,P) be the label of an edge in T\^{a,b). The (k — l)-th left 
child of the right child of this edge is also labelled by {a, (3). 

Proof. An elementary calculation shows that L = PDP^^, where 

fe'"/'' \ p _ A e"/'^' 

As a consequence, we get that for any integer j, 

1 / 




sin(7r/fcj \sin — ^ 
In particular, for j < k ~ 2, RU has nonncgativc entries, and 
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7r2 

7r4 




p^qV \p^q 



1 5 



2 4 4 2 8 

5 17 3 5 7 3 13 

Figure 2. The beginning of the tree R3(l, 1), which is a subtree of Ta3(1, 1). 

Therefore, for j < fc — 2, the j-th left child of the right child of the edge labelled 
by (a,/3) is labelled by {a,p)RU and the (fc — l)-th left child of the right child is 
labelledby (a, I -/3|). D 

We now define the tree Rfc(a, h) as the subtree of T^^. (a, b) obtained by removing 
the left child of the initial edge and all the edges which are (fc — l)-th left child. 

Proposition 2 (Linearity in Rfc(a, b)). Whenever it exists, the left child of an edge 
in Rfc(a,6) labelled by [a, (3) is labelled by {a,P)L. 

Proof. By definition, Rfc(a, 6) contains only j-th left children, for < j < fc — 2. 
The proposition is a direct consequence of the fact that for such j's, RL^ has 
nonnegative entries. D 

Considering Rfc(a, 6) as a subtree of T\i^{a,b), any edge in Rfc(a, 6) inherits its 
weight from T\^{a,b). As for T\{a,b), we organize the edges of Rfc(a, 6) in rows, 
and denote by 7r„ C ipn the set of edges in row n in Rfc(a, b) (see Figure [2|). 



In Section [5l we will first consider M(7r„), which is easier to study than 7\/ (?/;„). 
Then, the next step (Section [S]) will be to estimate M{'tpn) by partitioning the 
tree Tx^{a,b) in infinitely many copies of trees 'Rk{is+i,£s+2), where {is) is the 
sequence of numbers read along the leftmost branch of T;v^(a, b): ii := a, £2 := b 
and is+i := |Afe4 - 4-i| for s > 2. 

Remark 1. Let us explain why there is a mistake in the argument given in [5] 
(which deals with the particular case fc = 3 and p = 1/2). The average growth 
rate is proved to be equal to some explicit value for two linearly independent pairs 
of initial values for the random Fibonacci sequence which are: (31,52) = (Ij'/') 
and {l,ip~^), where (p = {1 + \/5)/2 is the golden ratio. It is then asserted that, 
since the vectors (1, (p) and (1, (fi^^) are linearly independent, any pair (a, b) can be 
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written as a linear combination of{l,ip) and {l,ip^^) to get the tree 'T\^{a,b) as a 
linear combination of the trees T\^{l,ip) and T\.^{1, ip^^); The conclusion follows 
that, since both of these latter trees show the same growth rate, the growth rate of a 
random Fibonacci sequence does not depend on the initial values a and b (with the 
obvious restriction that ab ^ 0). In fact, this way of proving the theorem is wrong, 
as it can be easily seen by writing the equality (1,1) = ip~'^{l,(p) + ip~^{l,(p~^). 
The tree deduced from this linear combination is a tree with root 1 and child of the 
root equal to 1, but it does not correspond to 'T\^{1^ 1) (since, for example, the left 
grandchild of the root is equal to 2ip^'^ instead ofO). 

5. Average M(7r„) of row n of Rfc(a, 6) 

For any X C T\i^{a,b), and < i < fc — 2, we define X^ as the subset of X 
made of all its elements which are i-th left children (X° thus corresponds to right 
children). With the convention ttj"^ := 7r2 and 7r| = for i < fc — 3, we get that 
each node of the tree Rfc(a, b) which has only one child is a (fc — 2)-tli left child. 

5.1. A recursive formula for M(7r„). The next lemma gives, for any < i < 
fc - 2, M{ttI^) as a function of (A/(7r^_ J)o<j<fc-2 and (M(7r^_2))o<j<fe-2- 

Lemma 2. For any integer n > 4 we have 

M(7r,"J = XpM{7r„^i) +pM {7:^-2) -pqMi7r';-l), 

M(7ri) = AgM(7r,",_i)-TOAfK_2), 

Af«) = AgM«l\)-q2^f«-|), 2<*<fc-2. 

Proof. Consider an edge e in 7r„ whose parent has weight w and is labelled by 
(a,/3). 

Case 1: Assume e € 7r° (e is a right child). The contribution of e to M{tt^) is 
{\f3 + a)wp = Xp.Pw+p.aw. Observe that when e runs over tt^, its parent runs over 
7r„_i and brings a contribution f3w to M(7r„_i). Moreover, when e's parent runs 
over 7rjj_]^, e's grandparent runs over 7r„_2 and has contribution aw/p to A/(7r„_2)- 
When e's parent runs over U^r;f7r*j_;^, e's grandparent runs over ^iZo'^n-2 ^^"^ ^^^ 
contribution aw/q to M {u'^'q Trl_2) = M(7r„_2) - M{tt'^zI)- This proves the first 
equality of Lemma [2l 

Case 2: Assume e S tt^ (e is a left child). Its contribution to M{ttJ^) is (A/3 — 
a)wq = Xq.Pw — pq.aw/p. Observe that e's parent is in tt^_i (thus is a right edge), 
and brings a contribution f3w to M{Tr^_^). Moreover e's grandparent is in 7r„_2, 
has weight w/p, and its contribution to Ai^(7r„_2) is aw/p. When e runs over tt^, 
its parent runs over 7r°_]^ and its grandparent runs over 7r„_2. 

Case 3: Assume e € tt^ for 2 < i < fc — 2 (e is a left child). The contribution 
of e to M(7r^) is (A/3 — a)wq = \q.(3w — q^ .aw/q. Observe that e's parent is in 
7r^l\ (thus is a left edge), and brings a contribution (iw to M(7r^~\). Moreover 
e's grandparent is in tt^I^j; tias weight w/q, and e is its only left grandchild in tt^. 
The contribution of e's grandparent to M{ttI/_2) is aw/q. When e runs over tt^, 
its parent runs over tt^^j^ and its grandparent runs over tt^jLj- This ends the proof 
of the lemma. D 

In the proof of the preceding lemma, we only used the structure of the tree 
Rfe(a, b) and the linear relation linking the labels of the edges in this tree, but not 
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the specific value A = Afc. In the next lemmas, this specific value plays a central 
role. 

Lemma 3. We have for any n > k + 2 

Proof. Consider an edge e in tt^Zi with weight w and label [a, /3). Its parent is in 
7r„Z2 and has weight w/q: Thus, it contributes for aw/q to Af(7r„Z2)- Moreover, 
e's (fc — l)-th parent is in 7r„_fc, has weight w/pq^""^ and its label {u,v) is such 
that {u,v)RL'^~'^ = (a,/3). Thus, (u,-y) = (/3, a — A/3) and this edge contributes 
for (a — Xf3)w/pq''~'^ to M(7r„„fc). The conclusion follows from the fact that when 
e runs over tt^Zi, its parent runs over tt^Zj ^-nd its (fc — l)-th parent runs over 

TTn-fc- □ 

Lemma 4. PKe /laue /or ani/ n> k + 2 

Proof. Consider an edge e in 7r„_fe with weight w and label (a,/3). If e G T^n-k^ 
j < k — 3, this edge is the ancestor of two edges in 7r^~^, which have labels 
{a, (3)RRL''~^ and {a, P)LRL''~'^, and weights wp'^q^~'^ and wqpq^~'^ . Their con- 
tribution to Af (tt^j"^) is thus (3wpq''~'^. If e S ^r^t-fc' ^^ ^^ ^^'^ ancestor of only one 
edge in tt^"^, having labels {a,P)RRL''~^ and weight wp^q^~^. Its contribution to 
M(7r^~^) is thus [3wp'^q^~'^ . The conclusion follows from the fact that any edge in 
7r^~^ has a unique ancestor in iin-k- D 

Lemma 5. We have for any n> k + 2 

A/(7r„) = AA/(^„_i) + (2p - l)M{Trn-2)+pq''-^M{7:r.-k) + {q - p)qM{7:'^zl). 
Proof. Using Lemma [2] and Lemma [31 we get 

fc-2 

Af(7r„) = 5^A/«) 

= (ApA/(7r„_i) + pAf (7r„_2) - pqM{Ti^-_l)) 
+ (AgA/(7r^i)-P'Z^/K-2)) 

A:-3 fe-4 

+ ^AQAf«„i)-^g2^.^(^«^_^) 
= AA/(^„_i) - \qM{TTl-_l) + {p-pq- <?')A/(^„_2) 

= AA/(^„_i) + (2p - l)A'/(7r„_2) + (q -p)qM {tt'^zD + pq'^-Hd {n,,^k) ■ 

D 

In the particular case p = q = 1/2, Lemma [S] gives M{'Kn) as a function of 
Af (7r„_i) and Af (7r„_fe). The next proposition gives a simple way to express Af (7r„) 
in terms of (Af (7rm)),„<„ in the general case. 
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Proposition 3. We have for any n > 2k + 2 

(3) M{Tr„) = AM(^„_i) + {2p - l)M(7r„_2) + pg*^"' AM(7r„_fc„i) 

+ p2q2fc-2M(7r„_2fc) 

Proof. Applying twice Lemma[5l (once for A/(7r„), tlien for M(7r„_fc)), the desired 
result follows from Lemma SI D 

5.2. Study of Af(7r„). 

Lemma 6. Let n > 2 and let Q{X) ~ X" — ^"jZq ^jX^ , where Cj > for all 
<.?<«. — 1. Then Q has a unique positive real root a. Moreover, a is of 
multiplicity 1 and all the other roots of Q have modulus strictly less than a. 

Proof. Since Q{0) = — cq < 0, Q has at least one positive real root a. Assume 
/3e*^ ^ a is another root of Q such that f3 > a. We have 

n — 1 n — 1 

j=0 3=0 

On the other hand. 



n-l 



^V' «.^....Y^ 



/3"-^ «" = E^./5- 



a J ^ — ' \ a 



n~j 



which gives a contradiction if /3 > a. Now, if /3 = a, we proved that 



^n-l 



E,=o Cjf3^e 



ijt 



S?=o Cj/9''. Since Cj > for all < j < ti — 1, this implies 9 = 2('!t and contradicts 
the fact that /?e** ^ a. 

If Q{ol) ~ Q'{a) = 0, then ?i X]?=o ^j'^'' ~ "Q^" = S?=o JCj'a-'i which is impos- 
sible since Cj > for all < j < n — 1. Therefore, a is a simple root. D 

Lemma 7. Let k > 3 be fixed and X ~ X^. The polynomial 

Pk{X) := X^^ - XX^^-^ + (1 - 2p)X'^''-^ - Xpq^-'^X^-^ -p^q^^'^ 

has a unique positive real root, denoted by ak, which is of multiplicity 1. 

Moreover, if p > Pc, all the other roots have modulus strictly less than a^. If 
P '^ Pc, Pk has 2 conjugate roots of modulus q, its positive root ai^ is smaller than 
q, and all the other roots have modulus strictly less than ak. 

Proof. We claim that Pk (X) can be rewritten as the product of X^ — qXX + q^ and 

fc-3 fc-2 

Y2k-2 ^„ Y2k-3 ^Y^„ „k-3~j yk-l+j ^2 V^ „ „2k-4:-j yj 

X - pa2k-3X -p 2_^ak-i+jq ■' X ■' - p 2_^ajq ■'A-' 

j=o j=o 

where all coefficients (aj)o<j<2fc-3 are positive. Indeed, identifying terms of degree 
j for < j < 2fc — 3 yields to the following relations, from which we can inductively 
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compute the coefficients 


af. 


Oo 


= 


1 


ai 


= 


A 


aj 


= 


Aflj-i — aj_2 V2 < J < k 


flfe-i 


= 


X+p{Xak-2 -flfe-a) 


ak 


= 


Aafc_i - pak-2 


aj 


= 


Aaj_i - aj_2 Vfc + 1 < j 


qa2k^3 


= 


Aa2fe-4 - a2fe-5- 



Observe that for < i < fc — 3, (a-;, a.i+i) = (1, \)U = (0, l)RU. Hence, we deduce 
from (HI) that a^ > for aU < i < fc - 2. Note that ^ also gives 

(4) i?i^-=(^ J 

This allows us to write {ak-z,ak-2)L ~ (0, 1)RL''~^ = (1,0), from which we get 
ak-2 = 1 and Xak-2 — o-k-z = 0. Hence ak-i = A > 0, and then afc = A^ — p > 0. 
Now, we have, for 1 < i < fc — 3, (a/j_2+i, flfe-i+O = {p,^)L'' = {qX,p)RL^. Thus, 
using ([2|) once again, qj > for all j < 2k — A. Next, using ([4]), we can write 

{a2k-5,a2k-i)L^ {qX,p)RL''^^ = (gA^+p, gA), 

and we get Aa2fe-4 — a2fe-5 = ^A. Hence we obtain a2fe-3 = A > 0, which proves 
that all the coefficients aj are positive. 

From the same equation above, we also get a2fe-4 = qX^ + p, and it is then a 
simple computation to check that, in the product, terms of respective degree 2A: — 2, 
2A: — 1 and 2fc also coincide with those of Pk- 

Since X^ — qXX + q^ has two conjugate (nonreal) roots of modulus q, we con- 
clude by using Lemma [HI that Pk has a unique real positive root ak , which is of 
multiplicity 1. Moreover, since Pk{q) = (7^'^~^(2 — A — 4p), we have as, > g if and 
only ifp> (2-A)/4. D 

Let us denote by ((3j) the roots of the polynomial Pk with ak ^ Poi and let /3i 
and /?2 = /3i be the two conjugate (nonreal) roots of X^ — qXkX + q^. Thanks to 
Proposition [3l for any n > 0, we have 

(5) Af (7r„+2) = Qo/3o" + E Qj(")/5r' 

where Qj is a polynomial (depending on a and b) of degree strictly less than the 
multiplicity of the root (3j . 

Moreover, since A/(7r„+2) is a real number, the coefficient Qq is real and coeffi- 
cients corresponding to pairwise conjugate roots are conjugate. 

Lemma 8. As soon as (a, 6) j^ (0,0), the coefficient Qq is positive. If p < Pc, the 
roots of X^ — qXkX + q^ are simple roots of Pk- If p < Pc, the associated coefficients 
Qi and Q2 are null. 

Proof. Assume ffist that p < pc- By Lemma[7l the roots oi Pk/iX^ — qXkX + q"^) are 
of modulus smaller than q. Hence /?i = ge** and /32 = ge~'^ (with 9^0 mod. 27r) 
are simple roots of Pk. We also know that the associated coefficients are conjugate 
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and that Po < q when p < pc- Hence, we deduce from ([5]) that when p < pc, Qi j^ 
impUes 

M(7r„+2) - 2q"gT:e(Qie™^). 

But this contradicts the fact that Af (7r„+2) is positive for all n. 

li p > Pc Po is the root of largest modulus of Pk- Moreover, since the coefficients 
Qi and Q2 are null whenever p < pc, we have M(7r„+2) ~ CqPq whenever p ^ pc. 
li p ~ Pc, we have Po = Qi thus 

M(^„+2) - q" (Qo + 25Re(Qi(n)e™^)) . 

If {a,b) 7^ (0,0), since M(7r„+2) is positive for all n > 1, we conclude that Qq > 
whatever p is. D 

6. Average M(V'„) of row ?i of Txi^{a,b), for A = Afe, A: > 3 

6.1. Left-most branch of Tx^{a,b). We denote by (^„)„ the sequence corre- 
sponding to the "left branch" of the tree T^^(a, 6), that is, the sequence defined 
as 

h:=a l2-=b in--=\Mn-i-tn-2\{n>i). 
Wc have the following result: 

Proposition 4. The sequence {£n)n is bounded. 

When Afc = 1 (fc = 3), it is easily seen that the sequence is upper-bounded by 
max(a,5). Strangely enough, the general case is quite more difficult to apprehend. 
Note that we cannot use Proposition [2] and L'^ = Id to say that this sequence is 
periodic (of period k) and thus bounded, since the ^„'s belong to the left branch of 
the full tree T\^. {a,b), which is not contained in Rfc(a, b). 

We give here a proof based on a geometrical interpretation, which can be applied 
for any < A < 2. This proof also appears in [5j. The key argument relies on the 
following observation: Let 6 be such that A = 2 cos 6*. Fix two points Pq,Pi on 
a circle centered at the origin O, such that the oriented angle {OPo,OPi) equals 
6. Let P2 be the image of Pi by the rotation of angle 9 and center O. Then the 
respective abscissae xq, xi and X2 of Pq, Pi and P2 satisfy X2 = Axi — xq. We can 
then geometrically interpret the sequence (£„) as the successive abscissae of points 
in the plane. 

Lemma 9 (Existence of the circle). Let 9 e]0, 7r[. For any choice of {x,x') £ 
R+ \ {(0, 0)}, their exist a unique R > and two points M and M' , with respective 
abscissae x and x' , lying on the circle with radius R centered at the origin, such 
that the oriented angle {OM,OM') equals 9. 

Proof. Assume that a; > 0. We have to show the existence of a unique R and a 
unique i g] — 7r/2, 7r/2[ (which represents the argument of M) such that 

Rcost — X and R cos{t + 9) = x' . 

This is equivalent to 

x' 
R cos t = X and cos 9 — tan t sin 9 = — , 

x 

which obviously has a unique solution since sm9 ^ 0. 

If X = 0, the unique solution is clearly R = x' / cos{9 — tt/2) and t = — 7r/2. 
Remark: Since xi > 0, we have t + 9 < 7r/2. D 
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Figure 3. i?„ = Rn+i is the radius of the largest circle, and R„ 
is the radius of the smallest. 



Proof of PropositionVR At step n. we interpret (.n+i in the following way: Applying 



the lemma with x 



and x' = In, we find a circle of radius i?„ > centered 



at the origin and two points M and M' on this circle with abscissae x and x'. 
Consider the image of M' by the rotation of angle 9 and center O. If its abscissa is 
nonnegative, it is equal to in+i, and we will have Rn+i = Rn- Otherwise, we have 
to apply also the symmetry with respect to the origin to get a point with abscissa 
H-n+i- The circle at step n + 1 may then have a different radius, but we now show 
that the radius always decreases (see Figure [3]). 

Indeed, denoting by a the argument of M', we have in the latter case 7r/2 — 9 < 
a < 7r/2, in = Rn cos a and £n+i = Rn cos{a + 9 + tt) > 0. At step n + l,we apply 
the lemma with x = Rn cos a and x' ~ Rn cos(a + 9 + n). From the proof of the 
lemma, if £„ = (i.e. if a = 7r/2), Rn+i = Rn cos(a + 9 + n)/ cos{9 — tt/2) = i?„. 
If In > 0, we have Rn+i ~ Rn cos a/ cos t, where t is given by 



cos 9 — tan t sin 9 



cos{a + 9 + n) 



cos a 



- (cos 9 — tan a sin t 



We deduce from the preceding formula that tant + tana = 2 cos 0/ sin 6* > 0, 
which implies t > —a. On the other hand, as noticed at the end of the proof of 
the preceding lemma, t + 9 < tt/2, hence t < a. Therefore, cos a < cost and 

Rn+l < Rn- 

Since £„ < Rn < Ri for all n, the proposition is proved. D 

Remark 2. The behaviour of the sequence (£„) strongly depends on the initial 
values a and b. R is proved in [5] that if a/b admits a finite X-continued fraction 
expansion, the sequence (in) is ultimately periodic. On the other hand, when £n 
decreases exponentially fast to 0, the exponent depends on the ratio a/b. Different 
examples of such behaviour are given in [2] . 
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6.2. A formula for A'/(-0„). For any s > and n > 2, we denote row n of the tree 
'Rk{£s+i,£s+2) by TTn^s- In particular, 7r„^o = T^n is row n of the tree Rfc(a, b). 

Proposition 5. For any n > 0, 

[n/fej n-km 

M{ll>n+2) = Y. Y^ CnMP<l'"'^rQ'M{TTn+2-s~km^s), 
m=0 s=0 

where c„,™ :== (J^) - (fc - !)(„" J. 

Proof. Let us code any trajectory from row 2 to row n + 2 in the tree T^^ (a, 6) by its 
successive steps {Xi)3<i<n+2 {^i = i? for a right step or L for a left step). Because 
of Proposition [U the label of the last edge of the trajectory does not change when 
we remove a pattern RL''~^ and its weight is divided by pq^~^ . By successively 
removing all patterns RL'^^^ in (Xi)3<i<„, we obtain a reduced sequence, that is 
a subsequence which never contains (fc — 1) successive L's, except eventually at 
the beginning. Moreover, its last edge has the same label as the last edge of the 
initial trajectory. Let us denote by s the number of left steps at the beginning of 
the reduced sequence and by m the number of removals. By removing these s left 
steps, we obtain a trajectory in the tree Rfc(£s+i,^s+2), whose last edge has the 
same label as the last edge initial trajectory and whose weight has been divided by 

Conversely, to a reduced sequence of length n — km (for < m < [n/fcj), we 
can associate a set of trajectories in T;^^ (a, b) ending in row n + 2 hy successively 
adding m patterns RL^~^ . Denote by Cn,m the cardinal of this set. To conclude the 
proof of the proposition, it remains to prove that Cn^m '■= (^) ^ {k — l)(m-i)- ^^ 
is obviously true for all n when m = 0. We claim that the sequence {cn,m) satisfies 
Cn+i,m = Cn,m + Cn,m-i for any n > and 1 < ?Ti < [{n + l)/k\. Indeed, there 
are two disjoint classes of trajectories in T\^. (a, b) of length n + 1 which can be 
obtained from a reduced sequence by successively adding m patterns RL^"^: 

• Those for which no pattern RL^^^ is inserted at the end. Their number is 
equal to Cn,m- Observe that this class is empty for n — km, since in this 
case, we start with an empty sequence. 

• Those for which at least one pattern RL''^^ is inserted at the end. They 
can be obtained by first inserting a pattern at the end, then by inserting 
(m — 1) other patterns anywhere but at the end of the new sequence. Their 
number is equal to Cn,m-i- 

It is straightforward to check that the numbers (J^) — (fc — l)(„,1i) satisfy the 
same induction, thus that they coincide with (c„^m). D 

Remark 3. The coefficients Cn^m can be obtained by generalizing Pascal's triangle 
(which corresponds to k ~ 1): We start by writing an infinite column of 1. On line 
k, we write a second 1 at the right of the first one. This new 1 is the beginning of a 
second column obtained by Pascal's rule : c„_|-i^„i — Cn^m + Cn^m-i- On line 2k, that 
is k lines after the beginning of the second column, we start a third column with the 
value C2k,2 = C2fc-i,i and go on using Pascal's rule. Each new column starts k rows 
after the previous one, its first term is given by the rule Cpk.p — Cpk~i,p-i (that is: 
the first term of a column is equal to the term one row and one column before it), 
and the next terms of the column are given by Pascal's rule. 
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Here is the beginning of the triangle in the case k = A: 

m = m ~ I m ~ 2 m = 'i m ~ A 
n = 1 

n = 1 1 

1 

1 

n = k 1 1 

1 2 

1 3 

1 4 

n = 2k 1 5 4 

1 6 9 

1 7 15 

1 8 22 

n = 3fc 1 9 30 22 

1 10 39 52 

1 11 49 91 

1 12 60 130 

n = 4fc 1 13 72 190 130 

1 14 85 162 320 



7. Proof of Theorem [T] 

Recall that we denoted by (/3j) the roots of the polynomial Pk (studied in 
Lemma [7]), with ak = (3o the eigenvalue with largest modulus. We deduce from ^ 
and Proposition [5] that there exists polynomials Qj.s (depending on is+i and is+2) 
of degree less than the multiplicity of f3j, such that for any n > 0, 

L"/feJ n-km 

j rn— s— 

We now study the contribution of each j3j : We will find an equivalent of the con- 
tribution of the dominant root f3o and prove that the contribution of other roots is 
negligeable. 

Lemma 10. For any j , the coefficients of the polynomials {Qj^s)s o-i"e uniformly 
bounded with respect to s. Moreover, the coefficients {Qo.s)s ore positive. 

Proof. This comes from Proposition [4] and Lemma [51 D 

\i p > Pc (respectively if p < Pc), let e > be small enough such that for all 
j > 1 (respectively j > 3), pj := |/3j|(l + e) < /3o- The previous lemma implies that 
there exists a constant K such that for all n, \Qj^s{'n)(3"\ < Kp^. The contribution 
of /3o to M{^n+2) can be written as t/„^o(/9o), where 

L^/^J / k—l \ ^^ n—km 

m=0 ^ ^ s=0 



GROWTH RATE OF RANDOM FIBONACCI SEQUENCES 



15 



i-— 1 \ ^^ n — km 

pq \ M—^ ( q 



The contribution of any other (3j can be bounded by KUn{pj), where 

Ln/fcJ 

m=0 ^ ■" ^ s=0 

Observe that UnS){x) is bounded by a constant times Un{x) since the coefficients 
{Qo,s)s are positive and bounded. 



E ^ 



f 



pq 



There exists C(x) > 0, 



Proposition 6. Let x > 0. We set f{x) := a; I 1 
depending only on x, such that 

• If X > q and x^ > {k ~ l)pq^~^ , then C/„.o(a;) ^ C{x)(f{x)) as n ^ oo; 

• Ifx>q then \Un{x)\ < C(x)(/(x))"; 

• Un{q) ~ 0{n) as n ^ oo; 

• If X < q, then Un{x) = 0(1) as n ^ cx3. 

Proof. Case 1: x > q. When j = (which corresponds to the dominant eigenvalue 
of Pfc), the coefficients (Qo,s)s are positive and bounded, so we can choose S large 
enough such that 

n — krn 

^X/ 



s=S+l 



<s. 



SetAix):=j:toQoA£ 



L(n-S)/fcJ 

A{x) x" E c„ 



Hence, 



pq 



fc-i 



m=0 

Observe that, since c„.„i 



Ln/fcJ 



(") - (fc- 1)( "J = (") "-M 



n\ n — fc??i+l 



pq 



k~l 



we can rewrite 



1 



pq 



fe-i 



E 



f - kAI/n + 1/n 
f — Al/n + 1/n '- -L ' JJ 



where A/ is a binomial random variable with parameters {n,9). Since M/n —* 9 
almost surely as n — > 00, we get that this expectation goes to -ji^ if 6* < 1/fc (that 
is, if x^ > {k— \)pq^~^), and the same is true if we replace n by n — S*. Therefore, 
if x^ > {k — 1)^(7'^^^, for n large enough. 



A{x)[\^-6\< 



x" 1 + 



Ur.fi{x) ( i-ke 



(1 + ^) 



pq 



fc-l 



< .t" 1 



pq 



fc-i 



For j 7^ 0, since x > q, \Un{x)\ is bounded above, up to a multiplicative constant 
C{x), by 

L"/feJ 

X / ^ Cn,r, 

because c„,m < ( "J . 

Case 2: x = q. We easily see that 

Ln/fcJ 

(7„(g) = Y. Cn.mP"?"-" (n - fcm + 1) < n. 

m=0 
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Case 3: x < q. 

which is, up to a multipUcative constant C(x), less than 

L"/fcJ 



D 



Proof of Theorem]^ Assume p > Pc- We know from Lemma [7] that f3o ~ ak > q, 
and all the other roots /3j are such that \/3j\ < ak- 

We use Proposition [SI and need to understand the variations of /. Elemen- 
tary computations show that f{x) decreases when x ranges from to Xmin := 
\/(fc — 1)^*^^^, then increases. Observe that f{q) = 1, hence /(a;,„in) < 1- 

We claim that ak > a^min- This is true if p < 1/fc because, in this case, q > .Tmin- 
This remains true if p > 1/fc: Otherwise, we would have q < ak < a;min, hence 
f{ak) < 1. Moreover, if |/3j| > q, we also have f{\(3j\) < 1, and the contributions 
of other /3j's is, by Proposition [6l at most linearly increasing with n. This would 
imply that Af (V'n) = 0{n). But we know from [2] that whenp > 1/fc, the n-th term 
of a (p, Afe)-random Fibonacci sequence almost surely grows exponentially fast. By 
Jensen's inequality, this is all the more true for its expected value M{ipn), so we 
get a contradiction. 

It follows from Proposition[n]that the contribution to M(-0„-|_2) oi (3o is Un,o{Po) ^ 
C(/?o)(./(/5o)) as 71 -^ oo, and that for j ^ 0, the contribution of /3j, which is 
bounded by KUn{pj), is negligeable with respect to UnfliPo)- This ends the proof 
of Theorem [T] in the case p > Pc- 

Assume that p = pc- Wc know from Lemma [7] that /?o = q, and all the other 
roots Pj are such that \(3j\ < q. Wc thus deduce from Proposition [6] that M('(/;„) 
grows at most linearly. 

Assume p < p^. By Lemma [51 we know that {Qi^s)s and {Q2.s)s a-re null. 
Moreover, we know from Lemma [7| that f3o < q and that f3j < (3q for all j > 3. 
Using Proposition [51 we conclude that M^ijjn) is bounded. D 

8. NON-ANALYTICITY IN THE NEIGHBOURHOOD OF 2 

Proof of CorollaryUi By Theorem[2l we know that for A > 2, t/(A) is a root of the 
polynomial Q\{X) := X^ — XX — (2p — 1). On the other hand. Theorem [1] says 

that for A = Afc and p> Pc, GiXk) = ak{p) 1 + ff(^ : where ak{p) is a positive 
root of Pk- Thus, for any fc > 3, we easily get that 

afe(p)''-' Qx, {g{\k)) = 2pq'^-^ {a\p) +re'-') > 0, 
which proves that, for any fc > 3, Q{\k) is not a root of Q\^. □ 

Proof of Corollaryl^ The growth rate of the expected value of a (1/2, A)-random 
Fibonacci sequence has a very simple expression: It is equal to A when A > 2 
(Theorem [2]) and to 2ak — Afc, where ak is the only positive root of the polynomial 
X'' — XkX''^^ — 1/2^^, when A = Afc. Indeed, in the case p = 1/2, the polynomial 
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PkiX) in Theorem [U can be rewritten as {X'' + l/2'')Qk{X), where Qk{X) := 
X^-\kX''-^-2~^. Moreover, Qk{ak) = imphes that a^ (l + ^) = 2afe-Afc. 

Observe that Qki^k) < 0, which obviously imphes that ak > A^. Since Qk{ctk) = 
0, we have al~^{ak-\k) = 1/2''. Thus, Uk > Afc > 1 proves that < Uk-Xk < '2~'^. 

Since Xk tends to Aoo = 2 when k goes to infinity, if ^'(2) exists, then we must 
have 

^'(2)^ Km fg(y-g(2)U 1+ lim 2^^^^. 

k^+oo y Xk — 2 J k^r+oo Xk — 2 

The numerator of the latter expression tends exponentially fast to 0, whereas 
the denominator is equivalent to 27r^/A;^, so we get G'{2) = 1. 

If Q is of class C^ at A = 2, let us write its Taylor expansion at order 2: 

g{Xk) = g{2) + (Afc - 2)g'{2) + i^i-^g"(2) + 0((Afc - 2f). 

We then get that g"{2) is equal to the limit of the ratio 2! 2(afc - Afc)/(Afe - 2)^, 
which is equal to 0. The nullity of the n-th derivative of g is obtained in a similar 
way, by an induction argument. 

Hence, provided g is of class C°° at A = 2, ^(2) = 2, g'{2) ^ 1 and 5(")(2) ^ 
for any n > 2. The only possibility for g to be analytic at A = 2 is to satisfy 
g{X) = A on a neighbourhood of 2. But we also have ^(Afc) = 2afc — Afc, which would 
imply ak = Afc for k large enough. This would contradict Pfc(Afe) = — 2^'^ < 0. D 

9. Open questions 

9.1. Critical value. Theorem [1] states that for p = pc, the growth of E(g„) is 
at most linear. The proof of this result uses the fact that the labels £„ on the 
leftmost branch of Ta^. (a, b) are bounded. It is proved in [2] that if a/b admits a 
finite A-continued fraction expansion, the sequence (In) is ultimately periodic. The 
arguments developed in the proof of Proposition [6] show that, in this case, E(5n) 
docs grow linearly. However, we believe that for most choices of the initial values 
a and b, these labels decrease exponentially fast to zero, and that this ensures that 
E(g„) is bounded. 

9.2. Numerical simulation. As suggested to us by Steven Finch, from INRIA, 
numerical evidence of the growth rate of the expected value of a random Fibonacci 
sequence is not easy to obtain. This is due to the different behaviour of g„ and 
of ]E(g„). For A = Afc, comparison with the result obtained in [2] shows that for 
^-^ < p < 1/k, the expected value of the n-th term of a random Fibonacci 
sequence increases exponentially fast, whereas the sequence contains almost-surely 
a bounded subsequence. When 1/k < p < 1, numerical estimation of the growth 
rate of the expected value of F„ given by Theorem [T] suggests that it is strictly 
greater than the almost-sure growth rate. This would imply that the variance of gn 
increases exponentially fast with growth rate at least twice the growth rate of the 
expected value. In [4j where the case p ~ 1/2 and A = 1 is considered, the growth 
rate of the variance is proved to be equal to 1 + -x/S. It would be of interest to know 
better about the exact value of the variance for any p and, more generally, about 
the moments of higher order. 

9.3. Generalization. 
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9.3.1. Linear case. It is also ol interest to consider the linear case, given by the 
relation gn+i = ^Qn ± 9n-i, where the signs are given by an i.i.d. sequence 
of Bernoulli random variables of parameter p. We then have the easy induction 
E(g„) = AE((7„_i) + (2j5— 1)E(5„_2) for any A. The corresponding polynomial has 
two real roots (A ± ^JX^ + 4(2p — l))/2, and we easily get an explicit expression 
of E(g„) depending on the initial values. The question of interest in this setting 
would be to study the exponential growth of E(|g„|). The almost-sure growth rate 
of |g„| in the linear case is studied in [2] for A = Afe (fc > 3) and A > 2, and turns 
out to be more difficult than in the non- linear case. The analysis of E(|g„|) is also 
more intricate. Although we can embed the tree Rfc(a, h) in the tree of all possible 
sequences, this embedding is more complex than what we describe in the present 
article: The main reason is that after the removal of a pattern RL^~^ , left and 
right children are exchanged. However we think our method can be adapted to the 
linear case. 

9.3.2. Underlying structure. In the present paper, the underlying probabilistic struc- 
ture is a Bernoulli scheme of parameter p. There is no significant doubt that our 
method extends without significantly new ideas to some more general processes, as 
for example the one in which the choice of the plus or minus sign is given by two 
coins alternatively tossed, the first one with parameter p and the other one with 
parameter p' . 

We may also investigate what happens with a deterministic rule, like the codage 
of an irrational rotation in the circle. This leads to some interesting constructions 
which involve substitutions. These are to be explained in a forthcoming paper. 
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